function deblurredImage = iterDeblurring(img, psf, iterations)
    % LUCY_RICHARDSON 使用 Lucy-Richardson 算法进行图像去模糊
    epsilon = 1e-7;
    estimate = 0.5 * ones(size(img)); % 初始估计图像
    psfMirror = rot90(psf, 2); % PSF镜像

    for i = 1:iterations
        estimateConv = conv2(estimate, psf, 'same');
        relativeBlur = img ./ max(estimateConv, epsilon);
        estimateUpdate = conv2(relativeBlur, psfMirror, 'same');
        estimate = estimate .* estimateUpdate;
    end
    deblurredImage = estimate;
end